(N 

O 

(N 



This line only printed with preprint option 

Incompressible Navier- Stokes Equation from Einstein-Maxwell and 

Gauss-Bonnet-Maxwell Theories 

Chao Niu and Yu Tian 

College of Physical Sciences, Graduate University of 
Chinese Academy of Sciences, Beijing 100049, China 

Xiao-Ning Wu 

Institute of Mathematics, Academy of Mathematics 

and System Science, CAS, Beijing 100190, China and 



\Q. Hua Loo-Keng Key Laboratory of Mathematics, CAS, Beijing 100190, China 



Yi Ling 

—^ . Center for Relativistic Astrophysics and High Energy Physics, 

r-| ! Department of Physics, Nanchang University, 330031, China and 



Institute of Mathematics, Academy of Mathematics 
and System Science, CAS, Beijing 100190, China 



> 

o 

■^ ■ Abstract 



r~^ ■ The dual fluid description for a general cutoff surface at radius r = Vc outside the horizon in the 

o: 

charged AdS black brane bulk space-time is investigated, first in the Einstein-Maxwell theory. Un- 
der the non-relativistic long-wavelength expansion with parameter e, the coupled Einstein-Maxwell 
^ ■ equations are solved up to ©(e^). The incompressible Navier-Stokes equation with external force 

h: 

5r ■ density is obtained as the constraint equation at the cutoff surface. For non-extremal black brane, 

the viscosity of the dual fluid is determined by the regularity of the metric fluctuation at the hori- 
zon, whose ratio to entropy density rj/s is independent of both the cutoff Vc and the black brane 
charge. Then, we extend our discussion to the Gauss-Bonnet-Maxwell case, where the incompress- 
ible Navier-Stokes equation with external force density is also obtained at a general cutoff surface. 
In this case, it turns out that the ratio rj/s is independent of the cutoff Tc but dependent on the 
charge density of the black brane. 



I. INTRODUCTION 



It has long been known that excitations of a black hole horizon dissipate like those of a 
fluid with viscosity r] = t^-q [li (see also (2|-|5|), which together with the Bekenstein-Hawking 



entropy density s = jj^ yields the dimensionless ratio - = j-- Interestingly, this fact has 



4tt' 



been related to similar results |6| in the AdS/CFT correspondence, which regards the change 



of radius r as equivalent to the renornialization group (RG) flow |7H15| and the case of black 
hole horizon corresponds to the IR limit of this flow. Especially, it is proved that the ratio 



does not run with the RG flow 
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25] . and so the universality of this ratio in both the 



horizon fluid and the standard AdS/CFT follows. 



According to the general holographic dictionary jl6|, the Brown- York tensor of the bulk 
gravity is dual to the expectation value of the stress-energy tensor of the boundary field 
theory Q. Under the lo„g-wavele,rgth hm.t. the boundary theory can be described by 
hydrodynamics, known as the gravity/fluid duality |18l. Il9l|. Especially, under certain non- 
relativistic limit (or scaling), the boundary hydrodynamics takes the form of the standard 



incompressible Navier-Stokes equation 



11 



27] . Recently, Bredberg et al give the precise 
definition of the boundary theory on an arbitrary cutoff surface r = Tc outside the horizon 
in the Rindler bulk space-time, which reduces the bulk gravitational dynamics to the incom- 
pressible Navier-Stokes equation on the boundary (20|. This framework is extended to the 
AdS black brane case 2l|, as well as in the Gauss-Bonnet gravity 22]. 



In this Letter, we consider an arbitrary cutoff surface r = Tc outside the horizon in 
the charged AdS black brane bulk space-times in both the Einstein-Maxwell and Gauss- 
Bonnet-Maxwell theories, further extending the framework of (20| and 22] . Under the 
non-relativistic long-wavelength expansion with parameter e, the coupled Einstein-Maxwell 
(or Gauss-Bonnet-Maxwell) equations are solved up to 0{e^). The incompressible Navier- 
Stokes equation with external force density is obtained as the constraint equation at the 
cutoff surface, with a viscosity rj satisfying the universality relation ^ = ^ independent of 
both the cutoff r^ and the black brane charge in the Einstein-Maxwell case. In the Gauss- 
Bonnet-Maxwell case, the incompressible Navier-Stokes equation with external force density 
is as well obtained at a general cutoff surface. However, it turns out that the regularity 
condition at the horizon determines the ratio 2 = ^{1 — 2{n — A)a[n — 1 — {n — 3)g^]} 
with a the Gauss-Bonnet coupling constant and n the space-time dimensionality, which is 



independent of the cutoff Vc but dependent on the charge density qh of the black brane. 

The rest of the Letter is organized as follows. In Sec. II, we present the metric and 
electromagnetic background and introduce the non-relativistic long-wavelength expansion, 
focusing on the Einstein-Maxwell case, where the bulk equations of motion are solved up to 
0{e^). In Sec. Ill, the dual fluid on the cutoff surface is analyzed with the incompressible 
Navier-Stokes equation obtained in the Einstein-Maxwell case. In Sec. IV, we extend the 
above discussion to the Gauss-Bonnet-Maxwell case. 

II. METRIC AND ELECTROMAGNETIC CONFIGURATIONS 

We consider the standard n-dimensional Einstein-Maxwell gravity with the action 

^ ^ d"xv^(i? - 2A) - i [ (Tx^/^F^^F^"", 



IGnG J ^ ^' '4 

where A = — ^"^ ^^ ' is the negative cosmological constant. The corresponding equations 
of motion are 

G^^ + A^/.^ + SttGT^^ = 0, 

V^F^^ = 0, (1) 

where T^;^ = ^g^,^Fp„FP" — F^pF^^P is the stress-energy tensor of the electromagnetic field. 
Our background is the charged black brane solution 

■2 2m g2 



ds^ = -f{r)dt^ + 2drdt + r^dx''dx'', f{r) 



72 j^n—3 ».2n— 6 ' 



A = J r\r. Adt (2) 

under Eddington-Finkelstein coordinates, where the index a runs from 2 to n — 1, ?Ti is 
the mass parameter and Q the charge parameter. The corresponding electromagnetic field 
strength is 

F = \ ^^- -—^dt A dr. 

V SttG r"-2 

For convenience, we take the AdS radius / = 1 hereafter. 

The induced metric on the cutoff surface r = r^ outside the horizon is 

dsl = -f{r^)dt^ + rldx^dx", (3) 



which is flat and kept fixed when perturbing the bulk metric. In order to introduce the fiuid 
degrees of freedom v"" (the velocity) and P (the pressure), two types of diffeomorphisms that 
keep (j3]) invariant are taken: 

1. Lorentz boost 

with ^'^ = -^^v"; 

2. Special rescaling 

r — )■ (1 — P)r, t — )■ 

of r, t and x"". 



^^''"^ rt, a;° ^ .^ \. x" (4) 



\m-P)n] ' (l-P)re 



Then we promote v"' and P to be the velocity field v"'{t,x) and the pressure field P{t,x), 
which makes the transformed bulk metric and electromagnetic field no longer be solution 
of the Einstein-Maxwell equations ([T]) . We also adopt the non-relativistic long- wavelength 
expansion parameterized by e — )■ and the scaling |20| 

dt ~ e^ 9„ ~ e, dr-r^l, P r^ e^ v" ~ e (5) 

with da = gfs:, under which the perturbed bulk Einstein-Maxwell equations can be solved 
order by order. Up to O(e^), the transformed bulk metric by both types of diffeomorphisms 



has been given in 



22| (actually for arbitrary /(r)) as^ 



r^fir) .„,„,. „ r:v" 



ds^ = -f{r)dt^ + 2drdt + r^dx^dx" - 2r\l - ' ^i^ { )v''dx''dt - 2^^dx''dr 

retire) J [re) 



^ /(^c) /(r-c) /(re) 

/('^) fire) f{rc) 

The bulk electromagnetic field should also be transformed by the above two types of 
diffeomorphisms. After promoting v"" and P to be (t, x)-dependent (but r-independent) fields 



^ Note that our /(r) corresponds to r'^f{r) in [22|, and we do not distinguish Va from w". 



and adopting the scaling (|5]), the perturbed electromagnetic field can be straightforwardly 
worked out as 



(7) 

up to O(e^). It should be noted that in this charged configuration, there are gravitational 

degrees of freedom (DoF) and electromagnetic DoF. If the electromagnetic DoF are taken 
into account, there is a boundary current dual to the bulk electromagnetic field (see e.g. 



25|). In our approach, we only consider the DoF induced by the above two kinds of (lifted) 
diffeomorphisms, which can be roughly regarded as gravitational, and do not turn on the 
independent electromagnetic DoF. That is, our focus is on the infiuence of the charged AdS 
black brane, as the background (compared to the uncharged one), on the properties of the 
dual fiuid, especially the Navier-Stokes equation and the viscosity t] (or the ratio -) appearing 
there. It turns out that the resulting dual charge density is constant at visible orders, and 
so the conservation law of the boundary current just coincides with the incompressibility 
condition of the dual fiuid (see (9) below) and does not give any new equation.^ 

Now we should substitute the perturbed metric ([6]) and electromagnetic field ([7]) into 
the Einstein-Maxwell equations ([T]) and see if they have already solved the equations up to 



0{e^). Different from the asymptotically flat case 23|, the perturbed metric (together with 
the perturbed electromagnetic fleld) only solves the Einstein(-Maxwell) equations ([1]) up to 
(9(e), even if we imposed the 0{e^) constraint equation daV"" = (incompressibility of the 
dual fluid on the cutoff surface). It is known in the case without electromagnetic fleld that 
a correction term 

r^J^{r){daV^ + dbv'')dx''dx^ (8) 



at O(e^) should be added to the metric ([6]) |22| . It turns out that this prescription also applies 



to our case, without the need of additional correction terms for the perturbed electromagnetic 
fleld (jTj). In fact, the Maxwell equations in ([1]) are automatically satisfled at 0{1) and 0{e), 
while at C(e^) lead to two equations 

dav'' = 0, :F'{r)dav'' = 0, (9) 



^ The independent electromagnetic DoF can be turned on, e.g. by lifting the black brane charge Q to a 
function of (t,x), which will result in a non-constant dual charge density at visible orders. 



which can be solved altogether by the incompressibility condition. Then from the Einstein 
equations in ([1]) at O(e^) there is essentially only one requirement 

r"-V(r)-F"(r) + r"-^[(n - 2)/(r) + r/'(r)]J^'(r) + (n - 2)r"-=^ = 
for -F(r) to satisfy,'^ which can be solved as 

Here the integration constant C can be determined as C = r}^~^ with r/j the horizon radius 
by the regularity condition of the perturbed metric at the horizon, provided our charged 
black brane is non-extremal and so /(r) only has a simple zero aX r = r^. The explicit form 
of J^{r) is irrelevant to the succeeding discussions, where the one more integration constant 
can be determined such that the induced metric ([3]) is kept fixed to all orders in e. 

III. DUAL FLUID ON THE CUTOFF SURFACE 

On an arbitrary cutoff surface r = r^ outside the horizon, there is firstly a thermodynamic 
description of the (equilibrium) fluid dual to the background configuration (^. In fact, the 



Br own- York tensor 



2^ 



^"^ ^ 8^^^^'^ ~ ^"^ ~ ^^"^^ ^^^^ 



on the cutoff surface, with Kij its extrinsic curvature and K = g^^Kij, is 

8nG r, ^f{~) 2 r, 

(11) 



25] . On the other hand. 



which is identified with the stress-energy tensor of the dual fluid 
the stress-energy tensor of a (relativistic) fluid in equilibrium is 

Uj = (p + p)uiUj + pQij (12) 

with p the energy density, p the pressure and m* = J- — (1, 0, ■ ■ ■ , 0) the normalized fluid 
four-velocity. Inclusion of the constant C in (ITO|) is equivalent to the replacement 

C C 



P^P- ^^, P^ P 



SttG' ' ' SttG' 



We have explicitly obtained the above form of equation for 4 < n < 10. 



which leaves the combination p + p invariant, so we can omit C if we only consider this 
combination.^ Comparing (ITT]) and f fT2|) . we find 



p + p 


1 

SnG^r 


,f'(rc) 
S 2 


re 


Noticing the entropy 


density 


Sc = 


1 rr' 




AG rr ^ 


on the cutoff surface, 


the local Hawking temperature 






Tr=- 


1 firn) 



IQnGJJ^) r' " 



(13) 



(14) 



the explicit form of /(r) in ([2]) and the horizon condition /(r/^) = 0, we have the familiar 
thermodynamic relation 

p + p- ScT^ = q^Pc (16) 

with (up to some unimportant constant factor) the charge density 

"« = ;^ (17) 

c 

and the corresponding chemical potential 

n-'2 . Q Q_. , . 

SnG^fir^) r^ r^ 

which is just a (red-shifted) electric potential difference. 

In the perturbed case ([6]), the Brown- York tensor can be worked out as 

SirGU.dx'dx^ = -^J{^)^-^^^^llI^de + ^|=(ZM + (^-3)/(re) ^^^,^^, _ ^^^2 



/(r-c 



c) 

'^)ydx''dt 






''' =[(n - 2)/(r,)P + ry]{l^y^de + -^^^(My^dx'^dx" 



2^/fW) ' ^' ' 2,/fi^)firc) r 



4 ^3 

]Pdx''dx" 



'' ■ '■' ■M)?-(n-l)M)--rMK 



2\/7W2/('-o) '■ 



2 ^c V " ^V ^2 ^c c-V ^2 ^cJ 






2v7W 



[1 + ]{re)r{re)\{dav' + dbindx^'dx" + 0{e'), (19) 



Note that C is essential for the regularity of tij as re —S- 00 26 1 



after imposing the incompressibility condition ([9]) in the 0{e^) part.^ The transverse com- 
ponents of the Einstein equations give the conservation law 

D%j = n^T^, = FjiJ\ r = -n^F^'' (20) 

of the Brown- York tensor, where D is the covariant derivative on the cutoff surface, n the 
unit normal of the surface, Fji the boundary electroniagnetic field and J* the boundary 



current dual to the bulk electromagnetic field (see e.g. 25|). In our case, the cutoff surface 
is always flat, so D is just d. For the perturbed metric ([6]), it can be explicitly checked that 



n^T^, = 0{e 



3^ 



So the leading order equation of the index j = t in fl20|l is 

dHu = ^^=(M)'^^^« = FtJ'^ = 

at O(e^), which is just the incompressibility condition ^. The leading order equations of 
the index j = a in 02Up are 



r 



4 



■^yjay"+y''aw'^+4ap)-^n- 



167rG/(r,)y7W ^ 

(21) 

at 0(e^), where we have defined 

^ - y^^^y^ - (- - 1)/(^^) - -^/(-^)(^)c/(^)c (22) 

and fa = FaiJ^ as the external force density (see e.g. (Sj). 

Now we can read off the viscosity from the Brown- York tensor ( TT9|) . In fact, the stress- 
energy tensor of a (relativistic) viscous fluid is 

tij = (p + p)uiUj + pQij - 2'qaij - C9{gij + UiUj) (23) 

with cTjj the shear and 6 = diU^ the expansion. We are only interested in the above stress- 
energy tensor up to 0{e^) under the non-relativistic expansion. To this order we have 6* = 
by the incompressibility condition ([9]), which renders the last term in fl23l) to vanish, and so 

^ij = -VidiUj + djUi), 



Our result above, with fire) presenting, is slightly different from (27) in [22|, where there is no f"{rc). But 
in the case considered in 
give the same result. 



22| (without electromagnetic field), we have checked that these two expressions 



where V means projection to the x° directions. Comparing (123|) with (fT9|) . we have 



1 + f{r,)r{r,) 
167rG 



1 r 



n-2 



IGttG rj: 



-2' 



which together with flT^ gives the cutoff-independent result 



T] 1 

Sn 47r 



Under the special rescaling (j3]) (with constant P) of the background configuration ([2]), one 
finds that the transformed energy density and pressure become 
{n-2)rl J{r)^, 



Ps = P + 



167rGJ/(r, 



,2 ^c 



P + 0{e' 



Ps = P 



rl f{r). 



/Wv ../Wv, 



167rG'y^/(re) 2/(rc) r 



{^y: -in- l)(^)'c - rc(^)'e']P + 0{e-^) 



respectively. So the ratio of pressure (or "pressure density") should be 22, |27| 



using flT3|) and fl22|) . Introducing the coordinates 



i=,/fW)t, F = r,a;^ 



(24) 



under which the induced metric on the cutoff surface becomes rjij, one can easily recast the 
conservation equation ( I2T]) as the standard incompressible Navier-Stokes equation 



~ 16nGJf{r,) 

dtP^ + /3U/3" + d^Pr - Vd^fi- = U fa = ^,;', fa 



r 



3( fir)y 
c V r'^ ) c 



(25) 



with external force density fa-, where the kinematic viscosity 

V = . 

P + P 



(26) 



IV. THE GAUSS-BONNET-MAXWELL CASE 



The action of the ra-dimensional Gauss-Bonnet-Maxwell gravity is 
I 



^ frxy^iR - 2A + aCan) - \ / rf"xv^P^,P'^^ 



^GB = R — ^R^yR^'' + R^yarR^"" 



with a the Gauss-Bonnet couphng constant. The corresponding equations of motion are 

Gf„y + Ag^,^ + aH^u + SvrGT^,, = 0, 

V^F^^ = 0, (27) 

where 



We again take the charged black brane solution 



28| 



2^ 



ds^ = -f{r)dt^ + 2drdt + r'^dx''dx'', 






/W = -(l-^/l-4a(l-:S^ + -Szi) (29) 



as our background, where we have defined a = [n — 3){n — 4)a for convenience. The 
electromagnetic field is of the same form as in the Einstein-Maxwell case. 

It turns out that the non-relativistic long-wavelength expansion for the perturbed con- 
figurations in the Einstein-Maxwell case, i.e. the perturbed metric ([6]) and electromagnetic 
field (jTj) up to O(e^), as well as the form (|8]) of the correction term to the metric ([6]), can also 
be used in the Gauss-Bonnet-Maxwell case. Again, the Maxwell equations are automatically 
satisfied at 0{1) and 0{e), while at 0{e^) lead to two equations 

daV'' = 0, J'\r)dav'' = 0, (30) 

which can be solved altogether by the incompressibility condition. Then from the Einstein 
equations at O(e^) there is again only one requirement for J^(r) to satisfy, which can be 
solved as^ 



r2(^^n-4 _ 2(n - 4)a[r"-5/(r-)]') ^ /(r) ' 
Here the integration constant C can be determined as C = t^~^ — 2(n — A)ar^~^ f'{rh) with 
r/i the horizon radius by the regularity condition of the perturbed metric at the horizon. 
The Brown- York tensor on the cutoff surface r = r^. in this case is 



22 



29 



301 



*u = ^:^^^^3ij - Kij - 2a{3Jij - Jgij) - Cgij] (31) 

with 

Jij = ]- i2KKikK^ + KkiK^'Kij - 2KikK'''Ki, - K^Kij), (32) 



^ We have verified this solution ioi 5 < n < 10. 
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where we have omitted terms that do not contribute for this flat cutoff surface with induced 
metric ([3]). It can be exphcitly worked out for the background metric (l29l) that 



iijdX dX r^ ^v 

SttG 



(n-2)f(r,),2&f(r,) 



/"■')"-^f^l=^ - ')'"' 



.^(/M , (!i^|/M _ .nr.C-t^^ . /|^) ,..V. 



-Cdsl . 
Comparing ( 133|) and (IT2|) . we find 



(33) 






The entropy density Sc and local Hawking temperature Tc on the cutoff surface are of the 
same form (IT4l) and ( ITSl) as in the Einstein-Maxwell case, which again leads to the thermo- 
dynamic relation flT6|) with charge density (fT7j) and chemical potential (ITS]) . 

For the perturbed metric ([6]), the Brown- York tensor on the cutoff surface can be worked 
out as 

(n-2)/(r,)^2a/(_re) 
3r2 



^TxGUjdx'dx^ = ^Jf{r,)^ ^ [ iy ^^)dt^ -^d^l 



rl ff'ir,) , (n-3)/(rJ . ^, .2(n-5)/(re) , f{r,)- 



+77=(^ . ^^i^^ - mr.i^^^^^ . ^) )..«..» 



'f{rc) ^ 



f2afir^ \ . 



1 - '-^) Un - 2)fir.)P + ry) -^M)'. 
J _ 2af{r,) \ rl v'v'fir) 



'dx' 
A / rl[rl^2af{r. 



^l ^2/7(^V2/(r,)(r2-2a/(r,)) 
{n - 1)(^); - r,(^)';)pdx'^rfx« 



\/2 



^2 ^c 



V ^ rl n-r r2 ^^V 2y7W 



+0(e^), (35) 
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after imposing the incompressibility condition (130|) in the 0{e^) part7 So the leading order 
equation of the index j = t in (120|) is 



dHu = (^^4^ - 1) "^^i^rAv" = n^T,, = 

at 0{e^), which is just the incompressibihty condition (l30l) . The leading order equations of 
the index j = a in (l20l) are 



-{ii-'-^){^yMv^+v%v^ ' 



d%, = "" , {{i-"-^^^^^){'-^)',{dtv'' + v'd,v'' + -daP) 

167rG/(r,)y7(^ 



/(re) ^^ _ 2a(^ + ^(^)'J ) [1 + /(rjr (r,)]9^^;n = /, (36) 



at O(e^), where we have defined 

and fa = n^T^a = FaiJ^ as the external force density. 

Now we can read off the viscosity from the Brown- York tensor ( l35ll . Comparing f l23ll with 
( 135|) . we have 

^ = ^^ f 1 - 2«(^^^ - ql)\ (38) 

with the charge density of the black brane q^ = -^=^, which together with 0141) gives the 
cutoff-independent result 

^ = i-(l-2a(^-gD) = ^(l-2(n-4)a[n-l-(n-3)g^]). (39) 



This ratio agrees with the known result for infinite boundary r^, — !■ oo |3lN34|. Under the 
special rescaling (jl]) (with constant P) of the background configuration ( p9|l . one finds that 
the transformed energy density and pressure become 



IQnGJfirr} 



2 



Ps = P 



r: 2a/(re) rg(rg + 2a/(r,)) /(r) ^ 



167rGy^/(r,) ^c 2/(re)(r2 - 2a/(re)) r 



-(^-l)(4?)'^-^c(7?)c]^ + ^(^')' 



^ Similar to the Einstein-Maxwell case, our result above is different from (57) in [23|, but we have checked 
that these two expressions give the same result when turning off the electromagnetic field. 
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respectively. So the ratio of pressure should be 



using ( !34|) and (137|) . Introducing the standard coordinates ( 124|) . one can again recast the 
conservation equation (136|) as the standard incompressible Navier-Stokes equation (l25l) with 
external force density. 

V. CONCLUDING REMARKS 

Under the non-relativistic long-wavelength expansion, we have solved up to second order 
of the expansion parameter the bulk equations of motion for Dirichlet-like boundary con- 
ditions at an arbitrary cutoff surface outside the horizon in the charged AdS black brane 
space-times in both the Einstein-Maxwell and Gauss-Bonnet-Maxwell theories, without con- 
sidering the independent electromagnetic DoF. The incompressible Navier-Stokes equation 
with external force density, as well as a cutoff-independent viscosity to entropy density ratio 
2, for the dual fluid on the cutoff surface has been obtained in both theories, while in the 
Gauss-Bonnet-Maxwell case the ratio | = ^{1 — 1{n — 4)a[n — \ — in — 3)g^J} depends on 
the charge density of the black brane. 

The framework used in our Letter seems to work well in all known cases of intrinsically 
fiat cutoff surfaces, but for general curved cutoff surfaces it can no l ong er be used. For 
some special curved cutoff surface, there have been some discussions J35|. Moreover, it is 
worthy to point out that very recently an alternative method has been introduced in 36 1 
by imposing Petrov type I condition on the cutoff surface. It turns out that imposing this 
boundary condition is equivalent to imposing regularity condition on the horizon at least in 
the near horizon limit such that the Navier-Stokes equation can be derived in a much simpler 
way. Following this approach the generalized framework applicable to the spatially curved 
space-time has been presented in (STf, and we expect that this strategy can be applied to 
our current work in future. 
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